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Rapeseed (Brassica napus L.) is grown in different geographical regions of the world. It 
is adapted to different environments by modification of flowering time and requirement 
for cold. A broad variation exists from very early-flowering spring-type to late-flowering 
winter cultivars which only flower after exposure to an extended cold period. B. napus 
is an allopolyploid species which resulted from the hybridization between B. rapa and 
B. oleracea. In Arabidopsis thaliana, the PEBP-domain genes FLOWERING LOCUS-T 
(FT) and TERMINAL FLOWER-1 (TFL1) are important integrators of different flowering 
pathways. Six FT and four TFL1 paralogs have been identified in B. napus. However, 
their role in flowering time control is unknown. We identified EMS mutants of the 
B. napus winter-type inbreed line Express 617. In total, 103 mutant alleles have been 
determined for BnCGFTb, BnC6FTa, and BnTFLl-2 paralogs. We chose three non-sense 
and 15 missense mutant lines (M3) which were grown in the greenhouse. Although only 
two out of 6 FT paralogs were mutated, 6 out of 8 BnC6FTb mutant lines flowered 
later as the control, whereas all five BnC6FTa mutant lines started flowering as the 
non-mutated parent. Mutations within the BnTFLl-2 paralog had no large effects on 
flowering time but on yield components. F1 hybrids between BnTFLl-2 mutants and 
non-mutated parents had increased seed number per pod and total seeds per plant 
suggesting that heterozygous mutations in a TFL1 paralog may impact heterosis in 
rapeseed. We demonstrate that single point-mutations in BnFT and BnTFL 7 paralogs have 
effects on flowering time despite the redundancy of the rapeseed genome. Moreover, our 
results suggest pleiotropic effects of BnTFLl paralogs beyond the regulation of flowering 
time. 
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INTRODUCTION 

Rapeseed (Brassica napus L.) is a major oil crop accounting for 
more than 60 million tons of seed and more than 20 million 
tons of extracted oil worldwide (http://www.worldoil.com/). This 
crop is widely cultivated in several temperate regions of the world 
such as northern Europe, Canada, China and Australia. Adapting 
flowering time to regional environmental conditions has been a 
major target of rapeseed breeding. A broad variation exists from 
very early-flowering spring-type to late-flowering winter cultivars 
that only flower after exposure to an extended cold period, a pro- 
cess known as vernalization (Iniguez-Luy and Federico, 2011). 
In rapeseed, flowering time and yield are closely linked to each 
other. Several genomic regions have been reported to contain 
major QTL for both traits. In a doubled haploid (DH) popula- 
tion derived from a cross between the Chinese semi-winter line 
Ningyou7 and the European winter-type Tapidor, at least four 
flowering time QTL were suggested as indicator QTL for yield 
(Long et al., 2007; Shi et al., 2009). Genetic variation within the 
different rapeseed types is relatively small, implying a need for 



wide crosses between non-adapted ecotypes to introduce traits 
of interest into elite materials (Girke et al., 2012). However, the 
introgression of genes from non-adapted ecotypes into elite cul- 
tivars is difficult due non-adapted flowering time. Therefore, it 
is of great interest to measure the effects of different paralogs on 
flowering and other characters. 

Brassica napus and its close relatives Brassica oleracea and 
Brassica rapa belong to the family Brassicaceae which also 
includes the model plant Arabidopsis thaliana. Rapeseed is 
an allotetraploid species originating from the spontaneous 
hybridization between B. rapa (AA) and B. oleracea (CC) less 
than 5000 years ago (Ziolkowski et al, 2006; Wang et al, 2011). 
These two progenitor species are ancient polyploids that under- 
went genome triplication between the Brassica- Arabidopsis split 
(~13 MYA) and their actual divergence event (~two MYA). 
Comparative mapping between the Arabidopsis and Brassica 
genomes revealed numerous homologous regions arranged in 
highly syntenic chromosome blocks. Many Arabidopsis genes are 
represented in the B. napus genome by at least three paralogs 



www.frontiersin.org 



June 2014 | Volume 5 | Article 282 | 1 



Guo et al. 



Mutations in single FT- and 7FM-paralogs of rapeseed 



(Schranz et al., 2006). Due to their close phylogenetic relation- 
ship and the high economic importance of rapeseed, knowledge 
transfer from the model species Arabidopsis to the complex 
Brassica genomes constitutes a worthwhile challenge for genomics 
research. 

In Arabidopsis, four pathways controlling flowering time have 
been extensively studied (Amasino, 2010). All these pathways 
converge at the CO/FT regulon (Andres and Coupland, 2012). 
Under long day (LD) conditions, the CONSTANS (CO) protein 
accumulates in leaves and induces expression of the floral integra- 
tor gene FLOWERING LOCUS T (FT) in the phloem companion 
cells (Moon et al, 2003; Turck et al, 2008). FT is the long-sought 
"florigen" and it is reported to be a strong mobile signal trigger- 
ing activation of floral identity genes in the Arabidopsis shoot 
apical meristem (Andres and Coupland, 2012). The FT protein 
is transported via the sieve tubes to the shoot apex, where it 
forms a heterodimer with the FD (FLOWERING LOCUS D) pro- 
tein (Abe et al., 2005; Wigge et al., 2005). Interestingly, a very 
closely related gene, TERMINAL FLOWER 1 (TFL1) plays an 
FT- antagonistic role by competing for FD, leading to a repres- 
sion of floral transition (Valverde, 2011; Andres and Coupland, 
2012). In Arabidopsis, FT-overexpressing plants and TFL1 non- 
sense mutants show the same early-flowering phenotype and 
produce terminal flowers in the shoot apex. TFL1 represses 
transcription of genes which are activated by FT (Hanano and 
Goto, 2011). In Arabidopsis, these two highly similar polypep- 
tides belong to a family of six members characterized by the 
phosphatidylethanolamine-binding domain (PEBP) (Kardailsky 
et al, 1999). Substitutions of crucial amino acids from the FT and 
TFL1 exon II, as well as the exchange of the exon IV led to con- 
trasting protein functions for both polypeptides (Hanzawa et al., 
2005; Ahn et al, 2006). 

Apart from their major role to control flowering time, FT and 
TFL1 orthologs have been shown to alter a variety of pheno- 
typic characters. In tomato, the SINGLE FLOWER TRUSS (FT 
ortholog) and SELF PRUNING (TFL1 ortholog) genes impact 
fruit yield heterosis. Fi hybrids generated by crosses between loss 
of function SFT mutants and tomato wild type ( WT) plants of the 
non-mutagenized line M82 have shown strong increment in fruit 
production (Molinero-Rosales et al, 2004; Krieger et al., 2010). 
FT/TFL1 gene orthologs have been characterized in diploid crops 
such as rice (Kojima et al, 2002), pea (Hecht et al, 2011), bar- 
ley (Faure et al., 2007), poplar (Bohlenius et al., 2006), and sugar 
beet (Pin et al., 2010) and in two polyploids, wheat (Yan et al., 
2006) and potato (Navarro et al, 2011). The characterization of 
the FT/TFL1 gene orthologs in polyploid plants is a special chal- 
lenge because duplicated genes can build new regulation networks 
leading to sub- or neo-functionalization (Pin and Nilsson, 2012). 

Rapeseed has six BnFT paralogs (BnA2FT, BnC2FT, BnC6FTa, 
BnC6FTb, BnA7FTa, and BnC7FTb) sharing high sequence sim- 
ilarity (92-99%) in their four exons (Wang et al, 2009). It has 
been shown that the BnC2FT copy is silenced in B. napus and 
B. oleracea due to the insertion of a miniature inverted-repeat 
transposable element (MITE) in its promoter region, whereas the 
remaining five copies are detectable in B. napus, B. rapa and B. 
oleracea (Wang et al., 2012a). There are at least four TFL1 par- 
alogs in the B. napus genome (Mimida et al., 1999). Among them, 



the BnTFLl-2 paralog shares high homology with the B. rapa 
ortholog on chromosome A10, whereas BnTFLl-1, BnTFLl-3 are 
highly similar to their B. oleracea counterparts. The BnC6FTa 
and BnC6FTb paralogs were co-located to a major flowering time 
QTL detected in nine winter-cropped environments which could 
support their function as flowering time regulators in B. napus 
(Qiu et al, 2006; Shi et al, 2009). Until now, B. napus FT/TFL1 
homologs have not been functionally characterized. 

This study had two major aims. First, we aimed to uncover the 
role in flowering time control of different FT and TFL1 paralogs 
in B. napus by analyzing EMS-treated offspring with missense 
and splice-site mutations within selected paralogs. We demon- 
strate that single mutations can change the onset of flowering 
in B. napus despite the redundancy of its allopolyploid genome. 
Moreover, we postulated that BnTFLl mutations also affect seed 
yield components in rapeseed. We found increased seed yield in 
Fi plants carrying a mutated BnTFLl allele on the Express 617 
background. Our data suggest that EMS-generated alleles may 
constitute a new resource to broaden the genetic basis of rapeseed 
breeding. 

MATERIALS AND METHODS 

PLANT MATERIALS AND GREENHOUSE EXPERIMENTS 

Seedlings of M3 lines and Express 617 (controls) were grown 
in the greenhouse at constant temperature (22° C) under long 
days (LD, 16 h light/8 h dark) for 4 weeks. Express 617 is an 
inbreed line (Fn) originated from the rapeseed winter-type cul- 
tivar Express (Harloff et al., 2012). Subsequently, plants were 
vernalized for 8 weeks at 4°C under LD conditions in a cold 
chamber. Of each M3 line, 30 plants were grown. After vernal- 
ization, plants were returned to the initial greenhouse conditions 
and transplanted to 11 x 11 cm pots. M3 plants and Express 617 
controls were arranged in randomized blocks. Plant positions 
on the greenhouse were indexed and linked to randomly gener- 
ated numbers using the Microsoft Excel software. Selected M3 
lines were crossed with the male sterile (MS) line MSL007 (NPZ, 
Hohenlieth, Germany) using homozygous M3 plants as a pollina- 
tors. Fi plants and Express 617 controls were grown in the green- 
house under the conditions mentioned above. F2 populations 
were produced by crossing M3 homozygous plants from a selected 
BnC6FTb mutation (BnC6FTbG2i54A) and non-mutagenized 
Express 617 plants. In each greenhouse experiment, the follow- 
ing phenotypic characters were measured according to the BBCH 
scale (http://www.jki.bund.de/en/startseite/veroeffentlichungen/ 
bbch-codes.html): first non-cotyledonal leaves (NCL, BBCH10), 
rosette plant (BBCH30), visible floral buds (BBCH50), first open 
flower (BBCH60), and end of flowering (BBCH69). Plants that 
did not grow beyond NCL (BBCH 10) were excluded from the 
experiment. Plant height, number of branches, initial flowers, 
filled pods, seed number and seed weight were recorded for each 
plant separately. 

MUTATION SCREENING 

A total of 3488 M 2 plants of the Express 617 EMS-population 
were screened by TILLING as described by Harloff et al. (2012). 
Gene specific primers were designed for BnC6FTa (FJ848915.1), 
BnC6FTb (FJ848917.1), and BnTFLl -2 (AB017526.1) 



Frontiers in Plant Science | Plant Genetics and Genomics 



June 2014 | Volume 5 | Article 282 | 2 



Guo et al. 



Mutations in single FT- and 7"F/_7-paralogs of rapeseed 



(Supplementary Table SI). For primer design and comparative 
analysis, B. rapa and B. oleracea genome sequences were down- 
loaded from (http://brassicadb.org/brad/downloadOverview. 
php) and (http://ocri-genomics.org/bolbase/), respectively. Plant 
genomic DNA arrayed in two dimensional 8-fold pools was 
amplified by direct or nested PCR. Forward and reverse primers 
were 5'-end labeled with 700 nm (DY-681) or 800 nm (DY-781) 
IRD fluorescence dyes, respectively (Biomers, Ulm, Germany, 
www.biomers.net). PCR amplifications with labeled oligos were 
done using the following profile: 95°C 5min; 35 cycles of 95° C 
30s, 60°C 45s, 72°C 90s, 72°C lOmin. Heteroduplex-specific 
restriction endonuclease CEL1 was extracted from celery and 
stored at — 80°C as reported by Frerichmann et al. (2013). 
Labeled fragments were separated by a LI-COR 4300 DNA ana- 
lyzer (LI-COR Biosciences) for 3:15 to 4:15h at 1,500V, 40mA 
and 40 W. Gel images were analyzed using the software GelBuddy 
(http://www.proweb.org/gelbuddy/). After positive pools had 
been identified, single plant DNA was amplified with unlabeled 
oligos and sequenced for SNPs confirmation. Sequences were 
analyzed with the CLC-bio Main Workbench sequence alignment 
tool (CLC bio, Aarhus, Denmark). 

DNA ISOLATION AND GENOTYPING 

Total DNA was extracted from young leaves using a CTAB pro- 
tocol (Morjane et al., 1994). Total DNA was treated with RNAse 
I (Fermentas, www.fermentas.de), and DNA concentration was 
determined by spectrometry (NanoDrop, www.nanodrop.com). 
DNA quality was checked by 1% agarose gel electrophoresis. For 
genotyping mutant lines, genomic DNA from single plants was 
amplified by PCR using unlabeled primers. PCR was done essen- 
tially as described in the previous paragraph. Five micro litter of 
each PCR product were loaded on 1% (w/v) agarose gels. Upon 
band size confirmation, the remaining 25 (xl of PCR product were 
sequenced via Sanger capillary sequencing. The sequences were 
analyzed with the CLC-Bio software (CLC bio, Aarhus, Denmark) 
using the sequence assembly viewer tool. 

TISSUES SAMPLING AND RT-qPCR 

Young leaves of M3 plants and Express 617 controls were sam- 
pled at four developmental stages, as described above. Genomic 
DNA sequences from the different flowering time genes ana- 
lyzed were retrieved from the non-redundant NCBI nucleotide 
database (http://www.ncbi.nlm.nih.gov/). Individual sequences 
were loaded to the CLC-bio main workbench version 6.0 (http:// 
www.clcbio.com), and groups of paralogs were aligned with 
the help of the internal alignment routine. Two main strate- 
gies were applied for expression analysis: (i) primers were 
designed in conserved regions within groups of paralogs for 
detecting joint gene expression, and (ii) copy-specific primers 
were designed for the members of selected paralog genes 
(Supplementary Table SI). Total RNA was extracted using the 
RNeasy kit (QIAGEN, www.qiagen.com) according to the man- 
ufacturer's protocol. The RNA concentration was determined 
by spectrometry (Nano Drop; Thermo Scientific, Wilmington, 
USA) and quality was checked by agarose gel electrophore- 
sis. Total RNA was treated with DNAse I (Fermentas Inc., 
Maryland, USA). First-strand cDNA was synthesized using 



Ohgo(dT)i8 primers and the M-MuLV Reverse Transcriptase 
(Fermentas). 

Quantitative real-time RT-PCR (RT-qPCR) was performed 
with SYBR qPCR Super mix w/ROX (Invitrogen Corporation, 
Carlsbad, USA) using a CFX96 Real-Time System (Bio-Rad 
Laboratories GmbH, Miinchen, Germany). Reactions were per- 
formed in a total volume of 15u,l containing 100 nM of each 
primer and 2 u,l of diluted cDNA templates, and amplified using 
the following cycling conditions: 95° C for 3 min, 40 cycles of 
95°C for 10 s, 60°C for 30 s, and 72°C for 30 s, followed by 95°C 
for 10 min. A melting curve was generated using a temperature 
range from 65°C to 95°C with increments of 0.5°C every 5 s. 
For each sample at least three technical replications were per- 
formed. For data analysis, the mean C t value of the target gene 
was normalized against the average Q value of two housekeeping 
genes (BnGADPH-3 and BnB-Tub). Calculation of relative expres- 
sion values was carried out following Pfaffl (2001) after extracting 
main Ct values via CFX manager software (Bio-Rad Laboratories 
GmbH, Miinchen, Germany). In each analysis, the relative expres- 
sion value for the reference sample has been set to 1. Normalized 
expression was averaged over two biological replicates and three 
technical repetitions in each case. Standard curves for the target 
and housekeeping genes are based on dilution series of purified 
cloned fragments for each gene. 

SEQUENCE DIVERSITY ANALYSIS 

For analyzing sequence diversity within the BnC6FTb and 
BnTFLl-2 genes, genomic DNA from one-hundred accessions of 
the B. napus ASSYST panel was amplified with paralog-specific 
primers and sequenced via Sanger method. We selected 117 lines 
from the B. napus ASSYST diversity set (Bus et al., 20 1 1 ) including 
winter, semi-winter and spring types which had been phenotyped 
in several environments worldwide (Supplementary Table S2). 
Lyophilized leaf samples harvested from young plants were used 
for DNA isolation with the NucleoSpin Plant II DNA isolation 
kit (Macherei & Nagel, Germany), following the manufacturer's 
instructions. PCR amplifications were carried out with paralog- 
specific primers as follows: 95°C 5 min; 35 cycles of 95°C 30 s, 
60°C 45s, 72°C 90s, 72°C 10 min. Sequences resulting from 
single band amplicons were assembled and aligned using the 
CLC-bio main workbench software (CLC bio, Aarhus, Denmark) 
and the resulting FASTA alignment was loaded into the soft- 
ware TASSEL (http://www.maizegenetics.net) for identification 
of polymorphic SNPs. 

RESULTS 

PARALOG-SPECIFIC EXPRESSION OF FIVE BnFT GENES 

We carried out a RT-qPCR experiment to measure the par- 
alog specific expression of six BnFT paralogs in leaves of 
the winter- type inbred line Express 617 during the transition 
to reproductive stages. Samples were taken from greenhouse- 
grown plants at three different stages of development (BBCH30 
before and after vernalization and BBCH50). Relative expres- 
sion values for each paralog were calculated after Ct normal- 
ization using BnGAPDH as a reference gene. Leaf samples at 
BBCH30 before vernalization (preV) were used as reference 
samples for relative expression calculations. At BBCH30 before 
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vernalization (BBCH30-preV), four BnFT paralogs (BnC6FTa/b 
and BnC6FT7a/b) were weakly expressed (Figure 1), whereas two 
transcripts were not expressed (BnC2FT and BnA2FT). Moreover, 
BnA2FT was only highly expressed at BBCH60 after floral tran- 
sition (data not shown), whereas BnC2FT showed no expres- 
sion at all. After vernalization (BBCH30-postV), BnC6FTa/b and 
BnA7FTa/b expression was higher in rosette plants, but differences 
between paralogs were obvious. BnC6FTb showed the largest rel- 
ative expression level (~9-fold). At BBCH50 (visible floral buds), 
the BnC6FTb and BnA7FTb paralogs showed the largest rela- 
tive expression levels (~13-fold) (Figure 1). In leaves at BBCH 
60 (first flower open), all paralogs with exception of BnC2FT 
showed very high relative expression levels (> 2000-fold), where 
BnC6FTa showed the highest relative expression (data not 
shown). 

EMS MUTATIONS IN BnFT AND BnTFLI PARALOGS 

We screened our EMS population to measure the flowering 
time effect of mutations within the BnFT paralogs BnC6FTa 
(FJ848915.1) and BnCSFTb (FJ848917.1). In Arabidopsis, apart 
from FT, other PEBP proteins such as TERMINAL FLOWER- 1 
(TFL1) regulate flowering by competing with FT for its bind- 
ing targets (Mimida et al., 2001). Therefore, we developed 
primers for the BnTFLl-2 (AB017526) gene assigned to B. 
rapa chromosome A10. In total, 3488 M2 plants were screened 
by TILLING for EMS-induced mutations in BnC6FTa/b and 
BnTFLI -2. 



We generated paralog-specific PCR amplicons covering 
between 50% (BnC6FTa) and 100% (BnC6FTb and BnTFLl-2) 
of the open reading frames. The BnC6FTa fragment covered 
exon I and intron I. The two BnC6FTb fragments covered 
exon I / intron I and exons III / IV (Figure 2). The BnTFLl- 
2 fragment covered all four exons. We identified 55, 14, and 
34 single nucleotide mutations in the BnC6FTb, BnC6FTa, and 
BnTFLl-2 genes, respectively. Forty-three mutations are located 
in introns, 19 are silent mutations, and three are located within 
the UTRs (Table 1). Mutation rates ranged between 1/72 kb 
and 1/24 kb per 1000 plants. The names of the mutant alle- 
les contain the nucleotide substitution and nucleotide position 
(Table 2). 

We identified one non-sense mutation in exon III of the 
BnC6FTb gene (5nC6Fn?Gi968A) leading to a truncated pro- 
tein by substitution of a tryptophan by a stop codon (position 
88). Another mutation (BnC6FTbc2009A) resulted in a splice- 
site deletion leading to a truncated protein by interrupting the 
junction between exons III and IV. The BnTFLI- 2 muta- 
tion in Exon III results in the substitution of a glutamine by 
a stop codon (position 146). Moreover, we detected numer- 
ous missense mutations in BnC6FTa (15), BnC6FTb (16), and 
BnTFLl-2 (10). 

We decided to focus on splice site- and missense-mutations 
which are most likely to affect the protein function. All observed 
missense mutations were compared to the SIFT database (http:// 
sift.jcvi.org) in order to evaluate the impact of the amino 
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FIGURE 1 | Relative expression of five BnFT paralogs in Express 617 
plants at three developmental stages before and after vernalization 
(dotted arrow). Plants at BBCH30 were analyzed before vernalization (preV) 
and after vernalization (postV). The time point BBCH30-preV was taken as 
reference for calculation of relative expression in all target genes. Two 
biological replicates and three technical repetitions were analyzed for each 
time point. Error bars indicate the standard error of the mean of the relative 
expression values. Expression levels of target BnFT genes were normalized 
against ESnGAPDH total expression. Ct values of the paralogs BnA2FT and 



BnC2FT were below the detection threshold at BBCH30 and BBCH50. 
Significant differences (P < 0.05) are depicted by asterisks. Differences in 
relative expression between the BnC6FTa I BnCSFTb (*) and BnA7FTa I 
BnA7FTb paralogs (**) at each time point were tested via f-test. Lines at the 
base of the bars indicate the comparison pairs. All samples were taken 
between zeitgeber 11 h and 12 h in each developmental stage. Express 617 
plants reached BBCH30-preV ~30 days after sowing. BBCH30-postV was 
registered ~90 days after sowing. BBCH50 was registered ~107 days after 
sowing. 
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FIGURE 2 1 Graphical presentation of 18 point mutations in two 
FT [BnC6FTa, BnCBFTb) and one TFL1 (BnTFLl-2) paralog of the 
B. napus winter-type cultivar Express 617. The exon (blue boxes) - 
intron (black lines) structure of each gene is shown. Locations of 



SNPs are depicted by the gray bar below. Numbers refer to the 
position of the mutations from the START codon. Only 
STOP/splice-mutations (red) and selected missense mutations are 
shown. 



Table 1 | EMS mutations in three flowering time genes detected by 
TILLING of the winter-type inbred line Express 617. 



BnC6FTa BnC6FTb BnTFL1-2 



Number of paralogs in rapeseed 


6 


6 


4 


Total number of M2 plants screened 


3488 


3488 


2092 


Sequence screened by TILLING (bp) 


1023 


1767 


1292 


Nonsense mutations 


0 


1 


2 


Missense mutations 


6 


15 


10 


Splice site mutations 


0 


1 


1 


Total number of mutations 


14 


55 


34 


Mutations/kb 


72 


30 


24 


M3 families selected for phenotyping 


5 


9 


5 


Total number of mutations 




103 





The BnCSFTb paralog was screened by two different fragments, 1021 bp and 
746 bp in size. 



acid substitutions on the protein function (data not shown). 
According to this analysis, we selected 18 M3 lines for growth 
experiments in the greenhouse (5 BnC6FTa, 8 BnC6FTb, and 5 
BnTFLl -2 mutations) (Supplementary Table S3) (Figure 2). 



PHEN0TYPIC CHARACTERIZATION OF BnFT AND BnTFL1-2 MUTANTS 

First, we confirmed the genotype of each selected M2 plant by 
Sanger sequencing. Then, M3 lines were grown in the greenhouse 
together with non-mutagenized Express 617 plants as a control. 
The phenological development of BnC6FTa and BnC6FTb lines 
was clearly different. All five BnC6FTa mutants flowered as the 
control, whereas six out of eight BnC6FTb mutants flowered later 
(Figure 3). The C6FTbci96&A mutants (stop mutation) showed 
a flowering delay of ca. 18 days, while C6FTbc20Q9A splice-site 
mutants flowered 29 days later as the control. Interestingly, 40 
and 54% of the C6FTbci96SA and C6FTbc2009A M3 mutants, 
respectively, did not bolt at all (Figure 4). BnC6FTb missense 
mutants started flowering 7 days (BnC6FTbc2UiA) up to 26 days 
(BnC6FTbGnA) later as the control. Apart from flowering time, 
reduced fertility was also apparent, mostly in BnC6FTb M3 plants 
(Supplementary Figure 1). 

To evaluate the effect of background mutations on flowering 
time, we produced an F2 population by crossing BnC6FTbc2i5iA 
M3 plants with non-mutagenized Express 617. BnC6FTbc2i5iA 
M3 missense mutants gave higher hybrid seed yield as the 
stop mutants and they flowered 15 days later as the control. 
BnC6FTbc2i54A M3 plants showed a reduced number of ini- 
tial flowers in comparison to other M3 mutants, however most 



www.frontiersin.org 



June 2014 | Volume 5 | Article 282 | 5 



Guo et al. 



Mutations in single FT- and 7FM-paralogs of rapeseed 



flowers were fertile. A total of 26 F2 plants encompassing all 
three genotypic classes were grown in the greenhouse together 
with Express 617. In agreement with M3 observations, homozy- 
gous F2 mutants (ft ft) flowered 13 days later than F2 siblings 



Table 2 | Nucleotide position and amino acid changes due to EMS 
mutations in 18 missense/non-sense mutations in three B. napus 
flowering time regulators. 



Gene 


Mutation 


Exon 


Amino acid 


Mutant code 








ci 1 Kct i+i itinn 

3UU0II LU llv/l 1 




BnC6FTa 


G37A 


Exon I 


Gly13Arg 


C6FTa Q374 




G52A 


Exon I 


Val17Lle 


C6FTa G52 A 




C74T 


Exon I 


Ser25Leu 


C6FTaci4T 




G104A 


Exon I 


Arg35lys 


C6FTaQ mA 




G163A 


Exon I 


Glu55Lys 


C6FTaQt 634 


BnC6FTb 


G17A 


Exon I 


Arg6Lys 


C6FTb G i 7A 




G124A 


Exon I 


Asp42Asn 


C6FTb Gm A 




C666T 


CArG Box 


CArG Box 


C6FTbc66BT 




G1968A 


Exon III 


Trp88Stop 


C6FTb G i968A 




G2009A 


Intron III 


Splice site 


CeFTb&oodA 




G2122A 


Exon IV 


Arg112Lys 


C6FTb G 2U2A 




G2133A 


Exon IV 


Gly116Arg 


C6FTb C 2133T 




G2154A 


Exon IV 


Val123Met 


C6FTbs2K4A 


BnTFL1-2 


G52A 


Exon I 


Val18Lle 


TFL1-2 Gb2A 




C518T 


Exon II 


Pro83Ser 


TFL1-2CBKT 




G750A 


Exon III 


Gly105Arg 


TFL 1-2ffj 5 o A 




G851A 


Exon IV 


Val108Met 


TFL1-2 G85 ia 




C965T 


Exon IV 


Gln146Stop 


TFL 1-2 C 965T 



homozygous for the wild-type allele (FT FT) which did not 
show any significant differences in flowering time as compared 
to non-mutagenized Express 617 (Supplementary Figure 2). 

The stop mutation BnTFLl- 2 did not lead to a major 
delay in flowering time. In contrast, the missense mutants 
BhTFL1-2qs2A and BnTFLl-2c750A flowered ~10 days later than 
the control (Figure 3). Since, the stop mutation is close to the end 
of the BnTFLl-2 gene, a functional protein may still arise after 
translation. Furthermore, BnTFLl&soA mutants exhibited mod- 
ifications in plant architecture which gave us a reason to select 
them for crossing experiments. BnTFLlc750A mutants developed 
normally during the early growth phase until reaching BBCH50 
(visible floral buds). The internode elongation phase was much 
longer as compared to Express 617, as a consequence, mutant 
plants were not able to stand by themselves after BBCH50. In this 
M3 line, the flower development limited the continuous growth of 
the floral meristem, whereas side branches continued flowering. 

THE EMS MUTATIONS IN BnC6FTb AND BnTFLl-2 PARALOGS ARE 
LOCATED IN HIGHLY MONOMORPHIC REGIONS OF EXON III AND 
EXON IV 

To investigate the genetic structure of those BnFT and BnTFLl 
paralogs with paramount impact on flowering time, we ana- 
lyzed the sequence diversity of BnC6FTb and BnTFLl-2 in B. 
napus by sequencing their complete exons III and IV in 117 5. 
napus inbreed lines from different geographic origins and growth 
types. Sequences selected for analysis in each gene after qual- 
ity trimming are deposited in Supplementary Table S2. While 
BnC6FTb exon III turned out to be highly conserved, exon IV 
exhibited larger sequence diversity. Within 41 bp exon III of the 
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FIGURE 3 I Flowering time point of 18 BnC6FT- and BnTFLl-2 mutants 
grown in the greenhouse at constant temperature (22°C), and LD (16 h 
light) after vernalization (4°C, 16h light, 8 weeks). Days to flowering 
(BBCH 60) was measured in M3 plants homozygous for the EMS allele. The 



non-mutagenized donor line Express 617 was used as a control. The number 
of plants analyzed is written in brackets. Differences in flowering time 
between homozygous mutants and control plants were tested via f-test. 
Significant differences (P < 0.05) are depicted by asterisks. 
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Express 617 C6FTa G37A C6FTa G52A C6FTa C74T C6FTa G104A C6FTa 



G163A 




C6FTb G17A C6FTb G124A C6FTb C66 6T C6FTb G2122A C6FTb 



J G2133A 




C6FTb G2154A C6FTb G i 968A C6FTb G20 o9A 




TFL1-2 G52A TFL1-2 C518T TFL1-2 G750A TFL1-2 G851A TFL1-2 C96 5T 

FIGURE 4 | Phenotypes of 18 B. napus BnC6FTa/b and BnTFLl-2 EMS M3 lines. Photos were taken as the non-mutagenized Express 617 plants started 
flowering (top left). Plants were grown in the greenhouse at constant temperature (22°C), and LD (16 h light) after vernalization (4°C, 16 h light, 8 weeks). 



BnC6FTb gene, only a single polymorphism was found at position 
2004 which corresponds to an allele frequency of 1%. The EMS 
mutation G1968A resides within a sequence domain which is 
monomorphic among all accessions investigated. For BnC6FTb 
(exon IV), six polymorphic regions were found with minor 



allele frequencies of <5.0%. The EMS-generated alleles (posi- 
tions G2122, G2133, and G2154) are residing in monomorphic 
sequences (Supplementary Figure 3). 

In contrast to BnFT genes, a higher variability in BnTFLl- 
2 exon III than in exon IV was found. In exon III the SNP 
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showing the largest variation was a C insertion at position 731. 
The EMS mutation G750A which is also located within this 
exon, was located in a fully monomorphic domain. For BnTFLl- 
2 exon IV, only a T/C polymorphism was found at position 
1030. In conclusion, there is a high degree of sequence con- 
servation within the analyzed sequences. Our EMS treatment 
created novel sequence variations within these highly conserved 
regions (Supplementary Figure 4). FASTA-formatted sequences 
for each gene are deposited in Supplementary File 1. All sequences 
have been submitted to NCBI (www.ncbi.nlm.nih.gov) (acces- 
sion numbers KJ533546 - KJ533625 and KJ533626 - KJ533728). 

A BnCff/Tfi-SPLICE-SITE MUTATION IMPACTS THE EXPRESSION OF 
OTHER FLOWERING TIME GENES IN LEAVES 

We reasoned that a loss of function of the BnC6FTb paralog 
directly impacts the transcriptional activity of other major flower- 
ing time regulators downstream of BnFT. To test this hypothesis, 
we selected the BnC6FTbc2009A mutant because, first the G2009A 
SNP causes a splice-site mutation that leads to a truncated pro- 
tein, and second, C6FTbc2009A M3 plants are characterized by a 
marked flowering delay of about 29 days compared to the Express 
617 control (Figure 3). 

We chose BnAPl and BnSOCl as putative downstream tar- 
gets of BnFT genes based on our knowledge from Arabidopsis 
(Yoo et al., 2005; Kaufmann et al., 2010). We measured their 
expression in leaves. BnC6FTbc2009A M3 plants were grown 
in the greenhouse under constant temperature and LD condi- 
tions. For expression analysis, young leaves of three different 
plants were taken at stages BBCH30 (pre and post-vernalization), 
BBCH50, and BBCH60. Arabidopsis API and SOC1 sequences 
were BLASTed against B. oleracea and B. rapa. High homology hits 
were aligned and primers were designed from highly conserved 
regions. Subsequently, joint expression of all paralogs was mea- 
sured by RT-qPCR. Gene expression levels of BnAPl and BnSOCl 
were normalized using BnGAPDH and BnB-Tub genes. 

We detected altered transcriptional activities of BnAPl in 
BnC6FTbc2009A mutants when compared to Express 617 con- 
trol plants. Control Express 617 plants at BBCH30 (preV) were 
used as reference sample for relative expression calculations. At 
rosette stages (BBCH30), BnAPl expression in the M3 mutant- 
was higher than in control plants, while at BBCH50 and BBCH60 
relative expression levels were 10-40% lower (Figure 5). Before 
vernalization (BBCH30-preV), BnSOCl expression in mutants 
was 2.5-fold higher than in Express 617. After vernalization, we 
detected a reduction of BnSOCl in BnC6FTbc2009A mutants com- 
pared to the control. When the first flower opened, the difference 
in expression between BnC6FTbc2009A mutants and Express 617 
was at its maximum. The altered expression in leaves indicates 
that a single BnC6FTb mutation may affect other major flowering 
time regulators. We expect that both genes are expressed in the 
shoot apical meristem as well. 

PERFORMANCE OF F 1 HYBRIDS USING THE BnTFL1-2 MUTANTS AS 
PARENTS 

In tomato, mutations in FT and TFL1 orthologs accounted for 
fruit yield heterosis in Fi hybrids (Krieger et al., 2010). We 
made an initial experiment to address the question whether 



B. napus orthologs might have a similar function. For produc- 
ing Fi hybrids, we selected BnC6FTbQ2oo9A and BnTFLl-2cjsoA 
homozygous M3 mutants as pollinators due to their late flower- 
ing phenotype in conjunction with an altered inflorescence (lower 
number of fertile flowers, Supplementary Figure l).We crossed 
homozygous M3 plants with the male-sterile (MS) line MSL007. 
The MSL007 line (MSL-Express) is an isogenic line of Express that 
carries the male sterility Lembke (MSL) genie male sterility sys- 
tem (Basunanda et al., 2010). Thus, no Fi heterosis was expected, 
except effects due to EMS mutations. 

Fi hybrids were vernalized and grown in the greenhouse with 
the parental lines and Express 617. Both Fi hybrids showed 
differential effects. BnC6FTbc2009A hybrids showed no signifi- 
cant differences in seed number per plant and total seed weight 
in comparison to Express 617 as the best parent. In contrast, 
BnTFLl-2cjsQA Fi hybrids had significantly higher number of 
seeds/plant (20%) and total seed weight (40%) as compared to 
the best parent (Figure 6). Although the effects of background 
mutations cannot be ruled out, these results could indicate that 
BnTFLl mutations impact heterosis in B. napus. 

DISCUSSION 

In the present study, more than 100 EMS-alleles have been found 
for three flowering time genes of B. napus. Based on previ- 
ous reports, the average EMS-mutation frequencies are expected 
to be lower in diploid species (~ 1/380 kb) than in polyploids 
(~l/50kb) (Till et al, 2007; Wang et al, 2012c). In this study, 
mutation frequencies ranged between 1/24 kb and 1/72 kb. Using 
the same EMS-population, mutation frequencies ranged between 
1/12 and 1/22 kb for sinapine biosynthesis genes (Harloff et al., 
2012). Although the observed BnFT /BnTFLl mutation frequen- 
cies are slightly lower, our results are in the range expected for 
polyploid species. The present mutants will be a valuable resource 
to study flowering regulatory networks in polyploids and they can 
be introduced into breeding programs. 

Our aim was to provide data about the function of FT and 
TFL1 paralogs in rapeseed. We found that, despite the redun- 
dancy of mutations in a single gene, either non-sense or mis- 
sense mutations in the BnC6FTb gene resulted in a marked 
flowering delay. This supports our hypothesis that BnFT par- 
alogs contribute differently to flowering time regulation. A large 
plethora of Arabidopsis reports on FT loss-of-function mutants 
have established a robust correlation between FT mutations and 
flowering time delay in Arabidopsis (Andres and Coupland, 
2012). Contrasting with expectations based on Arabidopsis TFL1- 
phenotypes, BnTFLl -2 mutants showed a slight delay in flowering 
time. In Arabidopsis, a single amino acid change in TFL1 (tfll- 
lcfyW5Asp) kd to early flowering and limited the development of 
indeterminate inflorescence by promoting the formation of a ter- 
minal floral meristem (Bradley et al., 1997). On the other hand, 
Fi hybrids derived from crosses between BnTFLl -2 M3- and rape- 
seed MS lines showed increased seed yield compared to BnC6FTb 
Fi hybrids and Express 617 controls. Thus, although the role of 
BnTFLl-2 involving flowering time regulation is not likely to be 
conserved compared to its Arabidopsis ortholog, TFL1-2 appears 
to be involved in yield-related traits as reported for its tomato 
ortholog SP (Jiang et al., 2013b). As we mention in the following 
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FIGURE 5 | Joined expression analysis of two BnFT downstream 
targets in the BnC6FTb G2 oo9 mutant (open boxes) and Express 617 
(filled boxes). Four developmental stages were analyzed before and 
after vernalization (dotted arrows). Plants at BBCH30 were analyzed 
before vernalization (preV) and after vernalization (postV). Two biological 
replicates (M3 plants) and three technical repetitions were analyzed for 
each time point. Error bars: standard error of the mean for biological 



replicates. Ct of target genes were normalized against the BnGAPDH 
and BnB-Tub total expression. The time point BBCH30-preV in Express 
617 control plants was taken as reference sample for calculation of 
relative expression. Differences in relative expression were pairwise 
tested (control Vs M3 line) via f-test. Significant differences (P < 0.05) 
are depicted by asterisks. All samples were taken between zeitgeber 
11 h and 12h. 



sections, confirming this hypothesis is a must for new research 
approaches. 

During evolution, duplicated genes may undergo dosage 
adjustments (Papp et al, 2003; Conant and Wolfe, 2008), non- 
functionalization, or sub-/neo-functionalization (Force et al., 
1999). We wanted to know whether different BnFT and BnTFLl 
paralogs gained different function by studying their phenotypes 
and their transcriptional activities. Through digital gene expres- 
sion analyses, differential expression within early generations of 
re-synthesized- (F1-F4) and natural B. napus accessions has been 
reported (Birchler and Veitia, 2010). Three highly similar genes 
encoding endoplasmic reticulum-bound sn-glycerol-3-phosphate 
acyltransferase-4 (BnGPAT4-Cl, BnGPAT4-C2, and BnGPAT4- 
Al) showed different expression patterns and altered epigenetic 
features (Chen et al., 2011) which is in accordance with the 
assumption that in polyploids orthologous genes are frequently 
expressed in a non-additive manner (Jiang et al, 2013a). We have 



also observed marked differences in the expression of six BnFT 
paralogs in support of the non/sub-functionalisation hypothe- 
sis. In regard to their position within a major flowering time 
QTL (Wang et al., 2009), BnC6FTb paralogs seem to play the 
most important role as flowering time regulators in winter type 
B. napus. More evidence has been given by the expression anal- 
ysis of putative FT downstream targets BnAPl and BnSOCI in 
late-flowering BnC6FTbc2009A mutants. SOC1 encodes a MADS- 
box transcription factor, acting as a floral integrator (Lee and 
Lee, 2010). The gene SOC1 gene is expressed in the shoot apical 
meristem, and SOC1 mutations lead to late flowering phenotype 
(Borner et al., 2000). However, SOC1 is also expressed in veg- 
etative organs (leaves) (Hepworth et al., 2002). A reduction in 
SOC1 mRNA was detected in the meristem of late-flowering ft- 
7 (Trpl38Stop) Arabidopsis Ler background (Searle et al., 2006). 
In the future, the activities of these genes shall also be studied in 
the shoot apical meristem. 
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FIGURE 6 | Growth experiments with vernalized F1 hybrids after 
crossing two mutants with a male sterile isogenic line (MSL007). Each 
Ft was obtained from two different crossing experiments. Mutant parents 
(M 4 lines), and Express 617 were used as controls. Yield components 
were determined on single plants grown in the greenhouse. The number 



of plants analyzed is given in brackets. Differences in flowers, seed 
number and seed weight per plant against Express 617 were tested via 
f-test. Significant differences (P < 0.05) are depicted by asterisks. Growth 
conditions: 22°C, 16 h light, greenhouse. Vernalization: 4°C, 16 h light, 8 
weeks. 



The phenotypic studies presented here gave further support 
to our assumption that BnFT paralogs do not contribute equally 
to flowering time regulation. In sugar beet (Beta vulgaris), two 
FT paralogs (BvFTl and BvFT2) were reported to antagonistically 
regulate flowering time (Pin et al., 2010). Knockdown of the FT 
potato paralog StSP3D resulted in a late flowering phenotype, 



while knockdown of the second paralog StSP6A had no effect 
on flowering time but on tuberization (Navarro et al, 2011). 
In line with these findings, we observed differential effects of 
BnC6FTa and BnC6FTb mutations. The strong effect of BnC6FTb 
as a flowering time regulator in rapeseed has been confirmed. To 
which extend BnC6FT& and other BnFT paralogs are involved 
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in flowering time control needs to be investigated in the future 
(e.g., by using other TILLING mutants). As determined by our 
expression analyses, the strongest case of non-functionalization is 
shown by the lack of expression of the BnC2FT copy. This result 
is in full congruence with Wang et al. report (2012a) where this 
gene copy was neither expressed in B. napus nor in B. oleracea. A 
series of recent studies has demonstrated that beyond flowering 
time control, FT-like proteins act as mobile or cell-autonomous 
proteins that mediate other developmental processes, such as 
growth, plant architecture, and tuber formation (Carmona et al., 
2007; Kinoshita et al, 2011; Navarro et al, 2011). In contrast 
to our BnFT mutant results, a previous analysis of sinapine 
biosynthesis mutants from the same EMS population as in our 
study, phenotypic or physiological effects had been observed 
only in double mutants (Harloff et al, 2012, Harloff, personal 
communication) . 

Although yield heterosis is regarded as a quantitative trait, 
single genes can contribute to heterotic effects through overdom- 
inance, such as the Arabidopsis Erecta locus (Moore and Lukens, 
201 1). As a first example of single gene overdominance, the yield 
of sft-4537/± heterozygous tomato plants was increased by up to 
60% in comparison to their parents after crossing high yielding 
M82 inbred plants with low-yielding homozygous loss of func- 
tion mutants (sft-4537) (Krieger et al., 2010). In our work, Fi 
hybrids carrying a BnTFLl-2 mutant allele had a higher seed 
yield as the Express 617 parent. Our study delivers the first 
insights about potential TFL1 -related heterosis in B. napus. In 
the future, experimental data are needed to verify this hypoth- 
esis. We tested our initial hypothesis by crossing mutants with 
the non-mutated donor line. Although, the data point at BnTFLl 
as a major gene for heterosis in rapeseed the possible impact of 
background mutations must be analyzed by additional hybrid 
combinations. Moreover, we will sequence the BnTFLl-2 loci 
from rapeseed lines with high and low combining ability (Qian 
et al., 2007). If our preliminary greenhouse data will be confirmed 
by field experiments, this study will offer new perspectives for a 
hybrid breeding strategy which make use of BnTFLl-2 sequence 
variations. 
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